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We derive QCD sum rules from the nucleon two-point function in nuclear medium, calculating its 
specral function in chiral perturbation theory to one loop. Our calculation shows the inadequacy of 
the commonly used ansatz to represent the two-point function as the sum of a nucleon pole term and 
a spectral integral over the high energy continuum. We also point out that it is the energy variable 
■^j- ' and not its square that must be used in writing the dispersion integrals. The sum rules are obtained 

by equating terms linear in nucleon number density of the spectral and operator representations of 
the correlation function. We use them to predict the curent-nucleon coupling, in good agreement 
with the value known from the vacuum sum rules. They also predict a definite deviation of the 
four-quark condensate from its ground state saturation value. 
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£■«•«. ; I. INTRODUCTION 

The extension of the original QCD sum rules [1] to finite temperature and chemical potential has been the subject 
of many investigations [2,3]. These were first extended to finite temperature by Bochkarev and Shaposhnikov [4], 
who pointed out the importance of continuum contribution in the low energy region. Working with the two-point 
f^S 1 function of vector currents, they showed that the contribution of the tttt intermediate state, though negligible in 
vacuum compared to that of the p-pole, is important at finite temperature. 

As shown by Leutwyler and Smilga [5], chiral perturbation theory provides a reliable method to calculate the 
correlation functions. They took the two-point function of nucleon currents at finite temperature and considered all 



the one-loop Feynman graphs of this theory to calculate the nucleon pole parameters at finite temperature. Later 



Koike [6] worked out directly the absorptive parts of these graphs to construct their spectral representations. He then 
evaluated the QCD sum rules in agreement with known results. 

The sum rules at finite temperature were subsequently taken up in Refs. [7,8], where their spectral sides were 
k>( ] calculated from the set of all one-loop Feynman graphs with vertices from chiral perturbation theory. It was shown 
that the contribution of each graph can, in general, be separated into a pole term for the communicating single particle 
and a remainder, that is finite at the pole. Clearly to evaluate the changes in the pole paramters, the remainder terms 
are not necessary. But in the evaluation of the sum rules, where we require their contributions at large external 
momenta, both the pole term and the remainder are of comparable magnitude. 

This observation questions the widely used ansatz for the two-point function given by a pole term (with modified 
parameters) together with continuum contributions from high energy region only. Such an Ansatz ignores the equally 
important contribution from the contiuuum in the low energy region and is therefore incomplete. 

In this work we investigate the QCD sum rules derived from the two-point function of nucleon currents at finite 
nucleon chemical potential [9,10]. Our work concerns the spectral side of the sum rules. As at finite temperature, we 
work out the spectral representation in nuclear medium in chiral perturbation theory to one loop. The contribution of 
the Feynman graphs are separared into the nucleon pole term and a continuum contribution in the low energy region. 
Again we find the Borel transform of these two types of terms to be of comparable magnitude. 

Another point overlooked in the literature is the absence of any symmetry of the amplitudes in the energy variable 
in the present case. Unlike in pionic medium (at finite temperature), the two-point function in nuclear medium has 
no density dependent singularities from the 'crossed channel', as we have no antinucleons in the medium. Thus if q 
is the external four- momentum and we restrict to = (qo = E, q = 0), there is no absorptive part of the amplitudes 
for E < 0. Clearly the amplitudes have no symmetry under E — > — E. As a result the dispersion variable must be E 
and not E 2 . 
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We also address the issue of sensitivity of the sum rules. The in-medium sum rules, as they are generally written, 
are actually the corresponding vacuum sum rules, to which are added density dependent contributions. If the vacuun 
sum rule is stable in a region of the Borcl parameter and the density dependent 'corrections' are small, as is usually 
the case, it is likely that the stability of the in-medium sum rules do not depend sensitively on these corrections. 

To restore the sensitivity of our sum rules, we equate terms linear in the nucleon number density on the spectral and 
the operator sides of the two-point function. Our in-medium sum rules are thus independent of the vacuum sum rule, 
being built out of the 'correction terms' only. We recall here that in the vacuum sum rules the unit operator provides 
the dominant contribution on the operator side, to which higher dimension operators add small 'power corrections'. 
In the present sum rules, the contribution of the unit operator, being independent of density, is absent altogether. 

In Sec. II we discuss the ingredients needed to write our sum rules. Here we present the form of the spectral 
representation given by the real-time field theory in medium. Then we derive the interaction vertices from chiral 
perturbation theory. We also list here the operators and their Wilson coefficients, that we shall include in expanding 
the operator product in the two-point function. The spectral representations of different Fcynman graphs are worked 
out in Sec. III. The sum rules are written and evaluated in Sec. IV. In Sec.V we comment on different aspects of our 
sum rules. 



II. PRELIMINARIES 



A. Two point function 

Although we restrict the sum rules to symmetric nuclear matter at zero temperature, our procedure to derive these 
is completely general. Thus we start with the ensemble average of the time ordered product of two nucleon 'currents', 

U(q)=i J rf 4 a;e i<zx (T?7(a;)r?(0)}, (2.1) 

where for any operator O, (O) = Tr [pO] /Trp ,p = e~P( H ~^ N \ H being the Hamiltonian of the system, (3 the inverse 
of temperature and N the nucleon number operator with chemical potential p. The nucleon current r](x)o,i, with 
spin and isospin indices D, i respectively, is built out of three quark fields, so as to have the quantum numbers of the 
nucleon [11]. In obtaining the spectral representation, we need only its hadronic matrix elements, but to expand the 
operator product in (2.1) into local operators, we must know its specific structure. Of the different possibilities, wc 
take here the preferred one [12], which for proton (i = 1) is 

VDtl (x) = e abc {u aT {x)Cru\x)) { l5l ^d c {x)) D , (2.2) 

where C is the charge conjugation matrix and a, b, c are the color indices. 

It is well-known that in a medium amplitudes take a Lorentz covariant form, if we introduce the four- vector u^, the 
four- velocity of the medium. Thus our two point function admits a Dirac decomposition, 11(g) = A(q)l+B(q)$+C(q)i/i. 
In what follows we write the sum rules for q = 0, when it reduces to II(go = E,q = 0) = A(E)1 + D(E)jq. 

We shall use the real time version of the field theory in a medium, where a two point function assumes the form 
of a 2x2 matrix. But the dynamics is given essentially by a single analytic function, that is determined by the 
11-component itself. Thus if -Fii(g) is the 11-component of the matrix amplitude, the corresponding analytic function 
F(q) has the spectral representation [13] 

F{E ' ^-nj E'-E-iE'e ^ ' (2 ' 3) 

As the chemical potential for particles and antiparticlcs arc equal and opposite, its presence in Eq.(2.1) leads to a 
loss of symmetry of the amplitude under E — > —E. The simplest example to observe this asymmetry is provided by 
the components of the 2x2 matrix of the free nucleon propagator. Denoting the nucleon field by ip(x), we have for 
the 11-component [14], 



jS(p) n = J d*xe^(T${x)$<p)) u = (j/> + m) 



p 2 — m 2 + ie 



{n-(w p )6( Po ) + n+(uj p )6(-po)} 2ir5(p 2 - m 2 ) 



(2.4) 



where n T are the distribution functions for nucleons and antinucleons, n T (uj p ) = {e l3 ( UJpT ^ + l} , ui p = \J m 2 + p 2 , 
m being the nucleon mass. The asymmetry is particularly explicit in the limit /3 _1 — > 0, in which we shall work, when 
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n~(ujp) — > 6(fi — uj p ) and n + (w p ) — ■* 0. As we shall see in the next section, for the two-point function (2.1) at q = 0, 
this has the consequence that the analytic function (2.3) has no imaginary part on the negative real axis in the E 
plane. The range of integration therefore extends only over to oo. 

Note that the analytic function S{E) corresponding to the free propagator-matrix in medium is independent of n T , 
coinciding with the free propagator in vacuum. But in keeping with the absence of singularities for E < from loop 
diagrams, we replace it by (for q = 0) 

S{E)^- ^ + m ^ 1 . (2.5) 

v ; 2m E-m + ie K ' 



B. Chiral perturbation theory 

To calculate the two point function (2.1), it is most convenient to use the external field method, where the nucleon 
current r](x) is coupled to an external spinor field f(x). The well-known QCD Lagrangian £q CD is thus extended to 

£ = £qcd + h + vf. (2.6) 

Chiral perturbation theory [15] gives the effective Lagrangian for the pion nucleon system in terms of the pion 
matrix field U(x) G SU(2) and the Dirac spinor ip(x) for the nucleon [16]. In what follows we ignore isospin symmetry 
breaking. At leading order, the effective Lagrangian in presence of the external field f(x) is given by, 

£eff + 4>{i$ - m)ip + -gA^ilb^ + £f , (2.7) 

where is the well-known Lagrangian for pions, whose interactions we do not need in this work. The matrix field 
u is given by 1 , u 2 = U and gA is the axial vector constant (gA = 1.27). £/ is the part involving the external field 
f(x). Taking note of the transformation properties of different fields, we can write it to lowest order as 

£ f = Xf(wp L + v)ip R ) + c.c, (2.8) 

where tpR t L are the right- and left-handed components of defined as V'fl.L = 5(1 ±75)^ an d A is the current-nucleon 
coupling constant. 

In the explicit representation U — exp(i<f> a T a / F^), where (f> a are the hermitian pion fields (a = 1, 2, 3), r a the Pauli 
matrices and F n the pion decay constant (F^ — 93 MeV), we have 

C eff = ^ - m)V - ^-J W T a m a + \f (j> - ^£75^^) + c.c. (2.9) 

We need two more vertices. One is the four-nucleon coupling derived by Weinberg [17]. Written in terms of the 
four component Dirac spinor if>, this nonrelativistic Lagrangian has the form 

V = [^(1 + 7o)V>] 2 " ^ W + 7o)775^] 2 , (2.10) 

where Cs and Ct may be obtained from the s-wave scattering lengths as 2 

Cs = -( ' V, G / 1 



506MeV J \52MeV 



x No confusion should arise from using the same u to denote the pion matrix field (only in this subsection) and the u-quark 
field as well as the four- velocity it M . 

2 Here we take the singlet and triplet scattering lengths as —16.4 fm and 5.40 fm respectively. These are obtained from np 
and pp scattering experiments [18]. Theoretical calculations based on chiral Lagrangians also reproduce these values [19]. But 
a different value for the singlet scattering length, namely —23.7 fm also exists in the literature [20]. Our sum rules give less 
satisfactory results with this latter value. 
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The other vertex is that of 77 — > ipipip. The most general form of this coupling would give r\ as 

v = ^Tkwtv)^ + 6 l (^rv^)r J r a v} , (2.11) 

i 

where the sum runs over all the independent Dirac matrices T l . For the particular Feynman graph we shall be 
interested in, only one linear combination of the ten coupling constants dj and bi will appear in its evaluation. 

C. Operator product expansion 

As is well-known, not only the Lorentz scalar operators but also the vector and tensor operators can have non-zero 
ensemble averages due to the availability of the velocity four- vector of the medium. Examples are the quark current 
operator qj^q = uj^u + d-f^d and the (traceless) energy momentum tensor operators Qf and Oj?„ of quarks and 
gluons respectively given by 

= qilf^D v q - ^-g^qq , (2.12) 



®% - -GU G " C + \9^G c a0 G af3c , (2.13) 

with the symbols having their usual meaning, (to is the quark mass in the SU(2) symmetry limit.) But as we already 
mentioned, the unit operator, that contributes dominantly to the vacuum sum rule, will not contribute at all to our 
in-medium sum rules. 

Thus the contributing operators of lowest dimension in the expansion of the operator product of nucleon currents 
are qq, qi/iq (=q'q in the rest frame of matter, = (1,0,0,0)). Then there are the operators of dimension four, 
namely (a s /7r)G a ^G fn ' a , @ f > 9 = u^UvQ^'S . Here we retain only the operator 0^ , as the coefficients of the remaining 
operators are small, arising from three contractions of the the quark operators in the two-point function. Among 
dimension five operators, the coefficient of qa^G^q turns out to be zero [7]. The other operators bring small 
contributions to the sum rules [10] and are ignored. 

Of dimension six operators, we retain, as usual, only the quark operators that have no derivatives or gluon fields. In 
the case of vacuum sum rules, their vacuum expectation values may be related to the square of the vacuum expectation 
value of qq, using the so-called factorisation or vacuum saturation. A similar approximation can also be made for the 
ensemble average, but simple model estimates seem to suggest that it may not be as good as in the case of vacuum. At 
this point we follow Ref. [10] to use factorisation and then replace the scalar-scalar condensate by the parametrisation, 

(qq) 2 -+ (1 - f)(0\qq\0) 2 + f (qq? , (2.14) 

where / is a real patrameter. 

The coefficient of operator product expansion at short distance are evaluated conveniently in configuration space. 
Such a method has been described in the Appendix of Ref. [7] for the case at hand. One first performs a Wick expansion 
of the two-point quark composite operators. The coefficients of the above mentioned operators are obtained from the 
single and double contractions of the quark fields. We state here the results of such a calculation, 

(Tv(x)fj(0)) - 4-e {<*«) + 2 U + (uiu) 



3tt 4 x 6 



+ {^(uu) 2 + 2x.u(uu)(uiu)} , (2.15) 



where we factorize the four-quark operators and omit the piece (ui/lu) 2 , being quadratic in nucleon density. Also we 
shall ignore the small effects of anomalous dimensions of the operators and assume their renormalization scale to be 
1 GeV. 



4 



III. SPECTRAL REPRESENTATION 



Here we want to calculate the (density dependent part of the) two-point function in the low energy region to first 
order in n, the nucleon number density in symmetric nuclear matter at zero temperature, 

- A f d3 P a< \ 2 Pf 

where pp is the Fermi momentum related to the chemical potential p by /i 2 = m 2 + pp. To this end we draw all the 
Feynman graphs with one loop containing a nucleon line. In addition to the single nucleon line forming a loop, we 
include also the loop containing an additional pion line: to account for singularities with the lowest threshold. They 
are depicted in Fig. 1 along with the free propagator graph (a). 




FIG. 1. Free propagator and one loop Feynman graphs for the two point function in the low energy region. The wavy line 
represents the nucleon current, while the continuous and the dashed internal lines are for nucleon and pion propagation. 

Any nonsingular term in the contribution of graphs is of no interest to us and will be omitted. Thus the free 
amplitude of graph (a) is given by 

n(^)( )=A 2 ^) = -^ T -i(l + 7 o), (3.1) 

on making the replacement (2.5). 

To evaluate the self-energy diagram (b) we rewrite the four- nucleon interaction term given by Eq. (2.10) conveniently 

as 

C ^ = 4E CrAB^CD^A^B^C^D), (3.2) 
i=S,T 

where V s = (1 + 70), and, T T = (1 + 70)775- Then the two point function for this diagram is given by 

n(E) (b) = -X 2 S(E)aS(E). (3.3) 
Here the self-energy a is a constant given by 

° = -\YlCij ^{-zHsip^nr + ps&OuP}, (3.4) 

where the tr(ace) is over 7-matrices. On inserting the density dependent part of Su from Eq. (2.4) it gives 
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<r=l(C s -Or)n(l+>Yo), (3.5) 



reproducing the result obtained earlier [21]. Eq.(3.3) then gives 

8m [E -m (E - m 



n (^)(6) = - — (^5 - C T )n { ^— + , F [> (1 + 70 ). (3.6) 



It is simple to evaluate the diagrams (c) and (d). Since two of the nucleon fields are contracted at the vertex itself, 
we may do so in Eq.(2.11) for r\ itself and retain the n dependent term to get 

r/(x) -> -(a + Pjo)nip(x), 

where a and (3 depend linearly on a» and &j. Then we get 

Tl(E) {c)+{d) = -2\n(a + p 7o )S(E) 
\ 2 (n 



E-m 



(l + 7o), (=(a + p)/X. (3.7) 



Here the coupling constant ( is unknown. We shall attempt in Sec. IV to find it from our sum rules. 

We now consider the diagram (e) with the ttN two-particle intermediate state. It is the prototype of all the 
remaining diagrams, so we derive its integral representation in some detail. Its contribution is 

n (3)(e) = ~jpf^ ( 3 - 8 ) 

where 

r(g)n = i f d A xe^ x -D{x) lllb -S{x) lllb . (3.9) 

The standard way to find the imaginary part of such a loop diagram is to write it in momentum space and integrate 
out the time component of the loop momentum, from which the imaginary part may be read off. But here this 
integration is difficult due to the presence of 9(±p ) in the nucleon propagator (2.4). We may, however, integrate out 
the po variable in the propagator itself, getting 

7^ )ll= / J^^ p [{^ + rn)(l-n-)e- i ^ + ^-m)n+e i ^}e(x ) 

- {{i> + m)n- e - lp - x + {j> - m)(l - n+)e ipx } 0(-ar o )l • (3.10) 

The analogous expression for the pion propagator is 

1 f d 3 k r i 

-D(x) 11 = j {2n)32LJk [ {(1 + n)e~ ik - x + ne ik - x } 9(x ) + {ne~ ik x + (1 + n)e ik - x ) 6(-x Q )\ , (3.11) 

where n is the pion distribution function at temperature n — (e' 3a)fc — 1) _1 ,uju = \J m 2 + k 2 . Though n + (ui p ) 
and n(u)k) are zero for the medium we are interested in, we retain them at this stage for generality and symmetry. 

With these expressions for the propagators, we may carry out both the x° and x integrations in Eq.(3.9) and find 
the imaginary part from the energy denominators, 

imr(g)n = 7T tanh(/3(<7o - ^)/2)ImT(q), (3.12) 

where 

ImT{q) = - J 7^4^/ (2^4^ & ~ m ) tt 1 ~ n ~ + n ) S (l ~ P ~ k) + (n~ + n)S(q -p + k)}. (3.13) 

Here the first term corresponds to 77 — > irN and the second to rj + 7r — > N, giving rise respectively to the discontinuity 
across the unitary and the 'short' cuts [22]. For q = 0, these extend in the i^-plane over E > m + and < E < 
m — m^. Note the opposite sign before n~ in the two terms; it ensures ImT ~ 0(n) for small n. 
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The complete result for ImT\\ has also another piece given by 

/d 3 p f d 3 k 
(2^7 j<^^ k (t + m ){{l-n+ + n)5{ q+P + k) + in+ + n)5( q+P -k)}, 

which is non-zero for E < only, but has no term proportional to n~ . Clearly this observation does not depend on 
the structure of vertices in the loop diagrams. We thus have the general result that there is no spectral function to 
one loop for E < in a medium with nucleons only. 

The three momentum integration in Eq.(3.13) can be carried out immediately to get 

ImT(E) = ±^M, lo = (E 2 + m 2 -ml)/2E , (3.14) 
E 

where 

/M = ^ (7o^ - m)n-(ui), 

The +(-) sign in front in Eq.(3.14) correspond to the unitary (short) cut. Inserting this result in Eq.(2.3) we get the 
spectral representation for T, 

V{E) -J E'(E'-EY (3 - 15) 

where the integration is understood to run over the two cuts for E > 0, with appropriate sign for f(E). 
Let us illustrate a convenient way to treat the integration in Eq.(3.15). We write it explicitly as 

rm r dE'f(E') dE'f(E') 

T(E) = L +m , VW~E) - Jo VW^E) ■ (3 ' 16) 

The range of the second integral can be mapped on to that of the first by the inverse transformation E' — > (to 2 — 
m 2 )/E'. Noting that lo and hence f(w) is form invariant under this transformation, we have 

rw -Cf /(g) (Fb ( " 7) 

Once we are on the unitary cut, the kinematics is determined by the first (5-function in Eq.(3.13). With q = 0, it gives 



E - v 7 ™ 2 + P 2 + Vml+p 2 , (3.18) 

so that E becomes the total energy of the ttN system in its centre-of-mass frame. The 3-momentum can be solved to 
give 

p 2 = {E 2 - (m + m^) 2 }-^ 2 - (m - to^) 2 }^ 2 . (3.19) 

The distribution function n~(ia) = 6(fi — u) restricts the upper limit of the integral to Ep given by Eq.(3.18) with 
p replaced by the Fermi momentum pp. Approximating f(E) to leading order in momentum, the 7-matrix structure 
in Eq.(3.17) factorizes and we finally get, 

T(E)=A(E)(l- 70 ), (3.20) 

with 

A{E) = -^i m+m , ^ {e^E (m* - ml)/E' - e) ' (3 ' 21) 

The spectral representation of the remaining diagrams can be obtained in the same way. For the vertex diagrams 
(f) and (g), we have 

n(«) (/)+(9) = 2X 2 S(q)A(q), (3.22) 
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where the imaginary part of A(q) is related to T(q) by 

ImA(q) n = -jj^ih + m)Mr(g)„. (3.23) 

TV 

We can again simplify such expressions by restricting to leading order in momentum. A further simplification results 
in the chiral limit (m, = 0), in which we shall evaluate the sum rules. With these approximations, the 7-matrix 
structure in A factorises as in T, getting 

A(E)=B(E)(1- 10 ), (3.24) 

with 

R(F) _ f dE'g(E') n-(u)3g A ( &-w?V 

B{E) ~]^-^> 9 W--8^4Fl{-^E—) ■ (3 - 25) 
With the replacement (2.5), Eq.(3.22) becomes 

K{E) u)+[g) = -^B{E){l- lQ ). (3.26) 

Finally the self-energy diagram (h) is given by 

n(«) (fc) = -\ 2 S(q)T,{q) S(q), (3.27) 
where the imaginary part of £(<?) is again related to that of T by 

ME(g) u = + m)/mr(g)n(rf + m). (3.28) 

It gives the integral representation 

£(£) = C(£) 7 o , (3.29) 

with 

, . f dE'h(E') , s n-<Lo)3g 2 A f E 2 -m 2 \ 3 
C(E)= h(E) = -^l-^ ——) . (3.30) 



E'-E ' y ' At: 2 AF 2 \ 2E 
Eq.(3.27) may now be expressed as 

So far the evaluation of graphs has yielded terms contributing to either the nucleon pole or the 7rA^ continuum. Now 
the last two terms in Eq.(3.31) appear as a product of the pole (simple and double) and the continuum integral. For 
interpretation and later convenience, we write them as a sum of such terms, getting 3 



t-t / 77i\ A 2 A 2 ( Cl /m + c 2 ci 1 f d 



dE'h(E') f dE'h(E' 



m){E'-E) J (E' — m) 2 {E' — E) 



(l + 7o), 
(3.32) 



where 



f dEh(E) _ 3g 2 A n f dEh(E) _ 

Cl ~ J E-m ~ 16F 2 ' C2 ~ J (E - mf ~ 



3 It may be obtained by expanding C(E) in a Taylor series with remainder about the nucleon pole. Alternatively, we may 
resolve the product of energy denominators into partial fractions. 
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This completes our evaluation of the Feynman graphs in terms of their spectral integrals. 

As a side result, we may now read off the modified nucleon pole parameters in nuclear medium. The dispersion 
integrals appearing above in the evaluation of graphs are all finite at E — to, though some are non-analytic, being 
proportional to fractional powers of n. Collecting the results for the simple and double poles, we see that the vacuum 
pole given by Eq.(3.1) is modified to 



A 



*2 



1 



£-m* 2 



o(l+7o), 



where 



A = A 



TO 



1 + 



3 
4to 



Cs — Ct + 



(3.33) 



At normal nuclear density, it gives to* — to = —322 McV, that is much too large compared to the value obtained from 
the self-consistent, relativistic Hartree-Fock calculation, namely to* — to = —87 MeV [23]. We postpone discussing 
this point to Sec.V. 



IV. SUM RULES 



To write the QCD sum rules, we need the Borel transform of the amplitudes with respect to a convenient variable. 
In the case of vacuum amplitudes, the dispersion variable q 2 in the space-like region is the the natural choice for it. 
With q fixed, we might think of replacing it by E 2 with E on the imaginary axis in the _E-plane. However, as we 
pointed out earlier, the loop graphs for nuclear medium do not generate any density-dependent singularity at all for 
negative values of E, precluding any symmetry of the amplitude under E — > — E. We thus choose Q = —E > as 
the variable with respect to which we take the Borel transform. It is then given by the differential operator, 



Bm = lim 

Q— >oo, n^oo 



1 



(n-iy. 



_d_ 



such that Q/n = M is held fixed, replacing Q by the new variable M . Thus 



B 



(Q 



i 



(k- 1)1 M k 



%jM 



It is now simple to find the Borel transforms of the contributions of individual diagrams evaluated in the last 
Section. To leading order in n we get, 

3 / 2rr?\ F 

^{n (fe) } = -(c s -c T )(i--)^(i + 70 ), 

F 

B M {tt( c )+(d)} = -C]^(l + 7o ) , 
3 F 

^M{n (e) } - ~7£p2 J^i 1 - TO) , 

B M {n w } = , (4.1) 

where we have introduced, for short, F — \ 2 ne~ m / M . The last equation shows clearly and simply the importance 
of the continuum contribution in the low energy region, that is neglected in the commonly used spectral ansatz for 
the two-point function. The self-energy graph (h) does contribute to the pole term shifting both the residue and the 
position as given by Eq.(3.33). As the Borel transform of the full contribution of this graph is zero, the continuum 
part must be equal and opposite to the pole term for large enough Q. 

We now turn to operator expansion. We specialize to the rest frame of the medium — (1,0) and for external 
momentum q^ = (— Q,0). Then the Fourier transform of the operator coefficients in Eq. (2.15) gives 
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n(Q) = -^Q 2 ln{-Q 2 ){{uu) +4(^)70) 



+ ^Qin(-Q 2 )(e f y f0 

+ ^((™) 2 7o + 2(uu)(u^u)). (4.2) 

The ensemble average of a local operator can be expanded in powers of nucleon number density. Thus for any 
operator R we have the familiar result to first order, 

d 3 p 



(R) = (0\R\0) + J ^-±^-^2{p,a\R\p,a}n-{po). 



where \p, a) is a one nucleon state of momentum p and spin and isospin denoted jointly by a and normalized as 

{p,a\p',a') = {2^f2 Po 8 a ^5\p-p'). 
If the constant, averaged matrix element is denoted by (p\R\p), it becomes 

(R) = (0\R\0) + 

We now apply this formula to the different operators appearing in Eq.(4.2). For the operator uu it gives 

(uu) = (0\uu\0) + £-n, (4.3) 

where the quark condensate in vacuum has the value, (0|mt|0) = (0|dd|0) = — (225MeV) 3 , and rh is the averaged 
quark mass m = (m u + m d )/2 = 7MeV. The so-called cr-tcrm 4 is a = m(p\uu + dd\p)/(2m) = AhMeV [24]. The 
matrix element for the operator v)u is also obtained, if we identify the quark current with the nucleon current as, 
u j nil + d-f^d = 3^7^ ip . We then get 

3 

(v)u) = -n. (4.4) 

Finally to obtain (6^), we need the nucleon matrix element of the energy-momentum tensor, 

(P|©£» = 2Af(p^ - \ gtlv p 2 ), 

where the constant A* is given by an integral over the nucleon structure function in the deep inelastic scattering. At 
Q 2 = IGeV 2 , it gives A* = .62 [25]. We thus get 

(0/) = '-mA^n. (4.5) 

Before we can take the Borel transform of Eq.(4.2), we have to pay attention to the singularities generated by the 
logarithms. Writing 

ln(-Q 2 ) = InQ + ln(-Q) = ln(-E) + InE, (4.6) 

we see that the first and the second terms give rise to the imaginary parts for E > and E < respectively. But 
we have no singularities for E < on the spectral side. So we must also remove the second term in Eq. (4.6). Then 
noting that 

B M {Q InQ} = M, B M {Q 2 lnQ} = -2M 2 , 

Eq.(4.2) immediately gives 



4 The <7-term in this Section cannot be confused with the self-energy a in Sec. Ill 
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B M {U(Q)} = ( JL + 3 7 o) + ^Mm 70 

+ ^(0|HO)(3 + /^ 7 o)}^ (4.7) 

The only other ingredient needed to write the sum rules is the continuum contribution at high energy on the 
spectral side. Following the usual procedure, we extract the imaginary part of the two-point function in this region 
from the logarithms present in Eq. (4.2) for the operator product expansion. The corresponding amplitudes are 
given by dispersion integrals starting from W, say. Being proportional to the operator matrix elements, these spectral 
contributions may conveniently be transferred to the operator side. 

It remains to equate the coefficients of (1 =F7o) m Eq s - (4-1) and (4.7) and include the continuum contribution from 
the high energy region to get the desired QCD sum rules, 

3A 2 e -m/M /MA 

16F 2 M 2 \ 9A m ~ 7 

and 

= 4 <0|a " |0 > ( f j, + 3 ) - ^ ( 3 + £) v * + h A ' MmV > ■ < 4 - 9 > 

where 

The two sum rules contain three unknown constants A, £ and /, besides W, the lower end of the high energy 
continuum. The coupling constant A is, of course, known from the vacuum sum rules [11], but we consider the 
possibility of its independent determination from our sum rules. 

As usual we look for an interval in the Borel parameter M, where the spectral and the operator sides of a sum rule 
overlap, leading to a plateau in the curve plotting the value of a 'constant' as a function of M. We vary W and / in 
the range .8 < W < 1.8 and < / < 1 respectively. Thus for each set of values of W and / in this range, we search 
for a reasonable plateau in the plots of A 2 and Q against M. 

We find that, in general, the constancy of A 2 extends over a wider region in M than for (. In this situation we first 
select values for the pair (W, f) for which there is such a constancy in A 2 . Next from among this set, we choose those 
for which the curve for ( also shows a reasonable plateau. The different curves so selected are shown in Figs. 2 (a) 
and (b). Our results for the ranges of different constants are 

W = (1.6 - 1.8) GeV, / = 0.2 -0.4, A 2 = (1.5 — 2.7) x M)~ 3 GeV 6 , ( = (20 — 65) GeV~ 3 (4.10) 

It is remarkable that our sum rule result for A 2 is so close to that obtained from the vacuum sum rules [11]. 
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V. DISCUSSION 



Instead of assuming an ansatz for the two-point function, as is usually done in the literature, we have calculated it 
in chiral perturbation theory. To one loop, the coupling constants appearing at the vertices in Feynman graphs are 
all known, except for £ arising from the vertex r\ — > (tp^tp . We determine it from the sum rules. 

As we saw at the end of Sec. Ill, the four-nucleon couplings lead to too big a mass shift (self-energy) of the nucleon in 
normal nuclear matter. It is related to the fact that the formulation of chiral perturbation theory for the two-nucleon 
system faces a problem, due to the proximity of the virtual or bound states to the theshold of NN scattering [17]. It 
makes the effective coupling constants rather large. Further, when multiplied by n for normal nuclear matter, these 
are still too large for the one-loop result to be valid [21]. A resummation of the loop diagrams with these vertices is 
supposedly needed [26]. 

We must, however, point out that this breakdown of the one-loop result in normal nuclear matter does not prevent 
us from using this four-nucleon coupling for writing our sum rules. The point here is that we are not working in 
normal nuclear matter. All we need is the validity of an expansion of the two-point function in powers of n for 
sufficiently small n. Since any resummation can only alter coefficients of terms quadratic and higher in n, our use of 
the coefficient of n to one loop remains valid. 

Our calculation does not support an assumption in the popular ansatz for the two-point function, namely that it 
can be well-approximated by the (dressed) nucleon pole term and an integral over the high energy region. We find 
that each Feynman graph, besides contributing possibly to the pole, also contributes an integral over the low energy 
region, which is not accounted for by the ansatz. Also there are no density dependent discontinuities of the amplituide 
for negative energy in nuclear medium. Hence the variable to write the dispersion integrals is E and not E 2 . 



12 



Finally we suggest a different use of the sum rules. There is a tendency in the literature to try to predict the 
nucleon self-energy from the sum rules. Indeed, this is in accordance with the works on the vacuum sum rules, where 
these are used to find the masses and the couplings of the single particles communicating with the currents in the 
two-point functions. However, it appears to us that the sum rules in the medium, as we have written, may be utilised 
more profitably to find the condensates in the medium. 

A way to carry out such a program would be to extend the sum rules to non-zero external three-momenta q. Here 
we retain the different four-quark operators as such, without invoking factorization. With the additional variable q 2 
at hand, it should be possible to extract numerical values for these oprator condensates in nuclear medium. These 
values may then be compared with model calculations [27]. 
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